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Abstract. The controversy concerning the temperature correction to the Casimir force has 
been ongoing for almost a decade with no view to a solution and has recently been extended to 
include semiconducting materials. We review some theoretical aspects of formal violations 
of Nernst's heat theorem in the context of Casimir Lifshitz thermodynamics and the role 
of the exponent of the leading term of the dielectric permittivity with respect to imaginary 
frequency. A general formalism for calculating the temperature corrections to free energy at 
low temperatures is developed for systems which do not exhibit such anomalies, and the low 
temperature behaviour of the free energy in a gap between half-spaces of poorly conducting 
materials modelled with a Drude type permittivity is calculated. 



1. Introduction 

The Casimir force [1], once merely a tiieoretical curiosity, is becoming the center of widespread 
attention in the wake of rapid developments in microtechnology. The enormous experimental 
progress made over the last decade towards accurately measuring this force [2-12] has created the 
need to calculate the Casimir force with high accuracy in realistic settings, taking into account 
such effects as material optical properties, surface roughness and geometry effects. Reviews of 
recent progress include [13,14]. 

It was realised quite early that an ambiguity existed as to the interpretation of Lifshitz' 
formula [15] for the Casimir attraction between dielectric half-spaces: when describing an ideal 
metal by taking the permittivity to infinity, different results were obtained at finite temperatures 
depending on the way the limit was taken. The ambiguity was originally sidestepped by 
prescription [16] and re-examined only much later by Bostrom and Sernelius [17] who concluded 
that due to the finite relaxation time of conduction electrons in a metal the transverse electric 
(TE) reflection coefficient of a metal-vacuum interface must vanish in the zero frequency limit, 
contrary to Casimir's ideal metal approximation in which reflection coefficients are set to unity 
at all energies. This was further supported by another study by H0ye et al. [18]. For finite 
temperatures this vanishing of the TE zero frequency reflection coefficient leads to a prediction 
of a relatively large reduction of the Casimir force between metal plates at finite temperature, up 



to 15% at 300K. The Bostrom-Sernelius analysis was opposed on thermodynamical grounds for 
violating Nernst's theorem (the third law of thermodynamics) which states that a nondegenerate 
system must have zero entropy at zero temperature [19]. Moreover it was concluded that the 
series of high accuracy experiments at Purdue exclude the thermal correction predicted by 
the theory in which the TE zero frequency mode does not contribute [10, 11]. The debate is 
summarised in recent reviews [20,21]. 

Recently an analogous ambiguity in the Lifshitz formalism was brought up for the case of 
semiconductors [22-25]. A formal violation of Nernst's theorem is once again the difficulty, 
this time due to discontinuous behaviour in the transverse magnetic (TM) reflection coefficient, 
whose value in the limit of zero frequency depends intimitely on the way the small density of 
conducting electrons in semiconducting materials at finite temperatures are taken into account. 

However, we do not expect such formal violations of Nernst's theorem stemming from the 
mathematical subtleties of the Lifshitz formula to have implications for the physics of the 
problem. In [18] it was concluded that on physical grounds, no TE zero mode should be 
present for real metals, and recently a quantum statistical mechanical treatment came to the 
same result [26]. For semiconductors, earlier statistical mechanical analyses by Jancovici and 
Samaj [27, 28] and by Buenzli and Martin [29] for ionic systems are of interest. What is 
found is that the effective separation between the plates increases as twice the ionic shielding 
length, which implies a non-local behaviour of the dielectric function. This increase in effective 
separation also means that the ionic contribution to reflectivity vanishes with vanishing ionic 
concentration. The results of refs. [27-29] are restricted to ionic systems, but in the high 
temperature (classical) limit they recover the ideal metal result corresponding to no TE zero 
mode. 

2. Formal violations of Nernst's heat theorem; general theory 

The Lifshitz formula expresses the Casimir free energy between parallel surfaces described by 
polarisation specific reflection coefficients where q G {j?, s} is the polarisation (assuming 
specular reflection and no coupling between p and s modes). For simplicity we shall assume the 
surfaces to be identical in the following, in which case the Casimir free energy at temperature 
T reads 



wherein a is the plate separation, p, s denotes TM and TE polarisations respectively, and \C,rn 
are the (imaginary) Matsubara frequencies so that Cm = 27rmT. As conventional, the prime on 
the summation mark signifies that the m = term be taken with half weight. We will be using 
natural units k-Q = h = c = 1 throughout. Henceforth we will frequently omit the subscript 
m on C,m and the various quantities depending on it. The integral in ([1]) is over all transverse 
momenta k_L of the field (_L denotes a direction parallel to the surfaces) and the substitution 
= + C,'^ has been made. 

In the case where the interfaces are between vacuum and a half-space made of dielectric 
material, the Fresnel reflection coefficients read 



Here e = e(iC) denotes the dielectric permittivity relative to vacuum. 

It is straightforward to verify that the values of in the limit C ~^ depend on the leading 
exponent of e(iC) as this limit is approached. For materials with mobile charges, models of the 
permittivity will typically diverge in the zero frequency limit, whereas that of a pure dielectric 
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isolator reaches a finite value, lim^^o£(iC) = Assuming e(iC) ~ iC/'^)^ {'^ is a constant) as 
C — > 0, one readily obtains the limits arrayed in table [T] in which 

l] <0. (3) 
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Table 1. Values of rq{iC 0,k) for 

_l Q I different exponents A. 

-2 fs(.K) 1 

< -2 -1 1 



The model permittivities which cause formal violation of Nernst's theorem have the common 
trait that the exponent A takes one value at all finite temperatures which changes abruptly at 
exactly T = 0. A general treatment demonstrates that such temperature dependence is necessary 
in order for a formal violation of the theorem to occur [30,31] as we will briefly explain. 

An example is the application of a Drude model to describe the dielectric response of an 
infinitely large and perfectly pure metal lattice, for which e(iC) is modelled as 

For a real metal sample of finite size, iy{T) reaches a nonzero value at zero temperature due to 
electron scattering on boundaries, impurities and imperfections, and A = — 1 for all temperatures. 
In this case entropy vanishes in the zero temperature limit as it should [32,33]. In a perfect lattice 
of infinite size, however, electron relaxation is solely due to scattering on thermal phonons, so 
that ~ as T — > 0. Thus A changes from — Ito— 2atT = 0, making jump discontinuously 
from zero to a finite value as seen in tabled! Clearly lD = Wp in ([3]) in the case A = — 2. 

Another example is the semiconductor whose conductivity vanishes as a function of T. If a 
Drude model is used to model the permittivity of such a material, 

a formal violation occurs when cr(T) vanishes at exactly T = 0. In this case A = —1 at all T 
until absolute zero, where it skips to A = and the magnitude of rp jumps discontinuously. 
Let 

lim To = Rn(K: A) (6) 

as tabulated (note that apart from the A = —2 s-mode, Rq is independent of k). It can be 
shown that when reflection coefficients jump discontinuously at C = 0, the free energy obtains 
a term linear in temperature equal to the difference of the m = terms of ([1]) as obtained with 
the two zero-frequency reflection coefficients respectively. If the leading exponent of e changes 
from Ai to A2 at exactly T = 0, therefore, it leads to a residual entropy =5^ = —d^/dT at zero 
temperature 



In particular, when Rq^K] Ai) = Rq{Xi) and Rq{K; X2) = Rq{X2) one may use the relation 



/ dKKln[l - R^e-^^"] = ^Li3(i?2) (8) 

Jo 



where Li„(x) is the nth order polylogarithmic function 
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Lin(-) = Ep (9) 
1=1 



to write 

^A,^A. = E {Li3[i?g'(Ai)] - Li3[i?2(A2)]} . (10) 

q=p,s 

In the particular cases of Drude modelled metals and semiconductors discussed above it 
follows immediately from tabled] and Eqs. ([7]) and pO|) that, respectively, 

1 

=5^_i^_2 = — / dKKln[l -f 2(^)6-2'"^]; (metals) (11a) 
4vr Jo 

= {C(3) - Li3[i?^(o)]} (semiconductors) (lib) 

where Rp{0) = (e — l)/(e + 1) and where we have used Li„(l) = C('^)) the Riemann zeta 
function. These two exponent transitions are those which come into play for metals and 
dielectrics respectively, or more precisely, upon plugging a Drude-type permittivity model with 
vanishing iy{T) or cr{T) into the Lifshitz formula and extrapolating to zero temperature. Other 
exponent transitions, naturally, would give other zero point entropy expressions. 
Note that by letting a) ^ 00 in ([3]) so that fs(K) —1, the entropy (lllaj) becomes 

^ (i-*00 C(3) . ^ 

which is the well known result for the so-called modified ideal metal model obtained by this 
procedure [18]. 



3. Free energy temperature correction for poor Drude conductor 

The model for the conductivity we will be studying in the following is assumed not to depend 
on T within a finite range of temperatures including T = 0, in which case it is clear from the 
above that Nernst's theorem will be satisfied. While the consideration of such a model cannot 
resolve such anomalies as reported in the previous section, it is nonetheless useful to establish 
benchmark results in various models within the Lifshitz formalism which exhibits very nontrivial 
behaviour in the joint limit of zero temperature and frequency C. 

We consider a semiconductor modelled by a Drude type permittivity such as ([5]), but where 
we assume a to be constant within a finite range of low temperatures including T = (the case 
a = was worked out in [22]). A more detailed treatment of this model may be found in [34]. 

The method used is to note that =^(0, T) depends on T only through the prefactor T and 
the Matsubara frequencies Cm = 27rmT, so we may write the free energy ([T]) on the forrr0 

00 

^ = r^/,(a)J]'9(M) (13) 

q=p,s m=0 
^ Note that this convention difTers sightly from that of [34]. 



where the function fq{a) is a convenient prefactor and we use the shorthand notation fi = mt 
where t = T/Tq is a dimensionless rescaled temperature to be defined in Eq. (llSp below. 

In the limit T = the sum becomes an integral, and we are interested in the difference 
between sum and integral, which may be found from the Euler-Maclaurin formula. It turns out 
g{n) is not analytical at m = so it is necessary to start the sum at m = 1 (or a higher value), 
writing 



r = 



V - / dm g{mt) = -^(0) - / g{mt)<lm + -g{t) - V -^g 
m=o Jo \ ^ Jo ^ fc=l^^''^' 

^5(0) - ^ g{mt)dm + ^g{t) - ^g'{t) + T^9"'{t) 



(2fc-l) 
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Since the correction terms in (|14p are evaluated at m = 1 and we are considering small T, we 
may choose Tq » T so that /i <C 1. We anticipate that when expanded for small fi, the function 
g{fj,) is of the form 



g{fi) ~ Co + ci/u + C3/X2 + C2i/i^ In(^) + C2/i^ + c^fi^ + /x ^ 



(15) 



Upon insertion into (|14p one finds that the cq and C2 terms do not contribute. Terms of integer 
powers of fj, are determined by a finite number of terms in the series ()14p . but for the terms 
containing logarithms or half- integer powers every term contributes. The series obtained is 
asymptotic but a meaningful value may nonetheless be assigned to all terms of (jlSp by defining 
the series by Borel summatiorH or zeta regular isation. The result for polarisation mode q = p,s 
is [34] 



A^, = r/,(a)r, 
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(16) 



giving terms proportional to respectively T'^,T2 , and T^. It is understood that the coefficients 
ci through C3 are polarisation mode specific. We shall content ourselves with expanding free 
energy to order in the present paper. 

Equation (|16p effectively reduces the problem to one of determining the expansion coefficients 
of Eq. psp . This task is still not trivial, however, and we consider only the "intermediate 
asymptotic" region in which conductivity is very small compared to inverse separation but 
much greater than temperature 

(17) 



T <C 47ra «; 



It is convenient now to define Tq = 2a, that is 



^ _ T _ Ci _ 27rT 
To 4:TTa Ana 



(18) 



Assumption (fT7|) ensures that f <C 1. The frequency which enters into is (i ^ T <^ a, so 
£{iC) simplifies to 

e = e+- (19) 
We will consider the p and s modes individually in the following. 



^ See appendix of [34] 



4. The TM mode 

The TM mode expression for Casimir Lifshitz free energy exhibits highly nontrivial behaviour 
near zero temperature. To simphfy matters we note from (llTh that the quantity 



a = 2a{ATT(T) (20) 

obeys a ^ 1. We will determine free energy corrections perturbatively in powers of a in order 
to obtain analytical results for the coefficients in (|15p . 
Substituting the integration variable 

x = 2na = (21) 

the TM free energy may be written 



{ATTcryt 



^ ( poo /A \3j. 

E'j/ dxxln(l-r^e-^)Uii^^,,(.n) (22) 

(.Jan J ^" " 



m=0 "^-^"/^ ^ ^ 

where gp{m) is now chosen to be the expression within the curly braces. To leading order in a 
one finds 

ln(l - r^e"^) = ln(l - ^^.e"^) + ff{a^); (23a) 



With the reflection squared coefficient now a constant with respect to x the integral in ()22p can 
be solved explicitly and expansion in /i as in (llSp yields the coefficients ci and C21 to leading 
order in a as simply [34] 

27r2 ^ ^ 

ci = — ; C2i = 8 (24) 

and with (jl6p the low temperature correction to free energy is to leading order in a in terms of 
T and a reads: 

TT^T^ C(3)r3 
72(47rcT)a2 ^{ATiaYa? 

An examination shows that the correction to this result is approximately a factor smaller as 
anticipated. 



5. The TE mode 

A similar procedure is performed for the s (or TE) mode. With the substitution 



X = . = ^; X = + (26) 



the TE term of the free energy ([T|) may be written 
^ (47ra)3t 



4^2 



^' dxxln(l-r,V"x^) =l^^'5.(m). (27) 



Again gs{m) is defined as the expression between the curly braces and in the new variables the 
squared refiection coefficient reads 



rl{x) = {yjx^ + l-xf. (28) 



Expanding the logarithm of (I27p to hnear order in a, 
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(29) 



allows us to determine the temperature corrections for these orders of a. Both terms on the right 
hand side of (f29l) give integrals over x which are explicitly solvable, and subsequent expansions 
in powers of fi give the coefficients 



ci 



--(21n2 - 1) 
4 
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a 
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(30) 



Clearly cs <C ci,C2z, being of linear order in a. The temperature corrections for the TE mode 

2 

thus read 



(47ra)r^ 
48 



r21n2 



1) + -V27rC(-|)a(47ra)2r 
b 



3- ~5 
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^^|^ + ^(ri). (31) 



This result is in fact in perfect agreement with that obtained for Drude metals [32] but includes 
one more order. Thus the concordance implies that a further expansion in a would only yield 
corrections of higher orders in temperatures, and that for the case of the TE mode the expansion 
in a was not essential for obtaining this result. Eq. ()3ip is therefore valid also when a is not 
small, as is the case for a good conductor. 

Notably the temperature corrections are independent of e to order for both polarisations 
and only enters in higher orders, an insight not included in [32]. 



6. Numerical verification of results 

The final results ()25p and (I3ip have been checked numerically to verify their correctness [34]. 
A plot of the theoretical correction ()25p compared to a direct numerical calculation of the free 
energy is provided in figured) 
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Figure 1. Free energy calculated 
using direct numerical calculation and the 
theoretical correction (|25p shifted to coincide 
with the numerical result at T = 0. 



Figure 2. Plot of the quantity R from Eq. 



A much more sensitive test is provided by the quantity 



R„ 



(32) 



where /S.J^^ is the theoretically predicted free energy correction ()25|) and A#'p "™ is that found 
by direct numerical calculation. The data used for calculation are (in SI units) a = lOOOnm, 
a^i/eo = 10^2s-\ e = 11.66 and loq = 8 ■ lO^^s'^. 

We have found that A^^^ is of the form A^*"" = -CT'^{1 - CiT) and assume A^^^ to have 
the expansion 

^^num ^ _DT'^{1 _ JJ^T + D2T'^ + •••), (33) 
which predicts the following expansion for R: 

R = - - D,)T -^[D, + C,{C, - D,)]T' + .... (34) 

In the special case where C = D and Ci = Di, this becomes 

R = -D2T^ + ff{T^). (35) 



Thus if the coefficient C is incorrect R would not converge to at T = 0, and an incorrect Ci 
would show as a linear behaviour at small temperatures. None of these effects are perceptible 
in the figure, demonstrating that the corrections to (j25p are small as predicted. 
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Figure 3. Free energy calculated 
using direct numerical calculation and the 
theoretical correction (I3ip to order shifted 
to coincide with the numerical result at T = 0. 



Figure 4. The difference between the graphs 
in figure [3] plotted against the term of ([3T]) . 



A similar high precision check of the result for the s mode (131 p was not possible with these 
numerical parameter values because the correction relative to the free energy at zero temperature 
is extremely small, of order 10~^ at IK. The correctness of the terms proportional to and 
T2 was however thoroughly verified for a good conductor in [32] . The term oc is numerically 
elusive because while requiring very high accuracy for verification when a <C 1, it is completely 
dominated by other terms for good conductors. A comparison of a direct numerical calculation 
with the prediction (I3ip to order in figure [3] reveals that the difference between these graphs, 
plotted in figure [U is in the same order of magnitude as the term in ()31|) , while the term 

5 

proportional to Ta is too small to be visible at this level. 



7. Conclusions 

We have reviewed the theory of formal violations of Nernst's heat theorem emphasising the 
way such a formal violation can only occur when the leading order behaviour of e(iC) with 
respect to ( undergoes a discontinuous change at exactly T = 0. Such apparent problems with 
the Lifshitz formalism occur when the double limit where T and ^ are both taken to zero is 
not unique and depends intimately on the exact way in which a material's dielectric (and, in 
general, magnetic) response is modelled in this limit. As a general remark, Nernst's theorem 
concerns zero temperature only, and it is not a priori clear that one can simply extrapolate 
between these two very different temperatures and use a result at one temperature to draw 
conclusions at another. In particular, if a system behaves essentially different at T = than 
at room temperature, a formal violation of Nernst's theorem by extrapolation is not necessarily 
worrisome. The results obtained complement those found in the case of zero conductivity in [22] 
and generalise TE mode calculations for Drude metals in [32]. 

Using a Drudc type model to describe a poor conductor whose conductivity stays finite at 
zero temperature we establish the low temperature corrections to Casimir Lifshitz free energy 
between two identical half-spaces separated by vacuum. As modelled, both TE and TM modes 
exhibit a quadratic temperature behaviour at low temperatures. 
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